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Abstract: We complete the survey for finite-source/point-lens (FSPL) giant-source events in 2016-2019 
KMTNet microlensing data. The 30 FSPL events show a clear gap in Einstein radius, 9 was < 0m < 26 uas, 
which is consistent with the gap in Einstein timescales near tg ~ 0.5 days found by Mrez et all in 
an independent sample of point-source/point-lens (PSPL) events. We demonstrate that the two surveys 
are consistent. We estimate that the 4 events below this gap are due to a power-law distribution of 
free-floating planet candidates (FFPs) dNprp/dlog M = (0.4 + 0.2) (M/38Mẹ) ? /star, with 0.9Sp<1.2. 
There are substantially more FFPs than known bound planets, implying that the bound planet power-law 
index y = 0.6 is likely shaped by the ejection process at least as much as by formation. The mass density 
per decade of FFPs in the Solar neighborhood is of the same order as that of ’Oumuamua-like objects. In 
particular, if we assume that "Oumuamua is part of the same process that ejected the FFPs to very wide or 


unbound orbits, the power-law index is p = 0.92+0.06. If the Solar System’s endowment of Tate 


objects in Neptune-like orbits is typical, which is consistent with the results of [Poleski et_al. 


(2021), then 


these could account for a substantial fraction of the FFPs in the Neptune-mass range. 


Key words: gravitational microlensing, asteroids 


1. INTRODUCTION 


The claim by of ~ 2 Jupiter-mass 
free-floating planet candidates (FFPs) per star opened 
up the field of FFP population studies. This led to 
several works that tried to interpret this population 
(Clanton & Gaudi, 2014alb,|2016). However, it was also 
later contradicted by (2017), who found no 
evidence for such planets in a larger sample of events. 
Both of these studies were limited to statistical analysis 
of the Einstein timescale (tg) distribution, for which the 
mass dependence is convolved with dependences on the 
lens-source relative proper motion (Hre1) and lens-source 


relative parallax (7.1). 
(2021) initiated a new approach to probing 
the FFP population that focused on analyzing the Ein- 


stein radius (g) distribution. This approach has two 
advantages. First, Hre, = 0g/tg is automatically deter- 
mined, which removes one of the two convolutions with 
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unknown distributions that was mentioned above. Sec- 
ond, the selection function is relatively independent of 
lens mass (M) over a broad range of masses. 


When reported KMT-2019-BLG- 
2073 as the fourth microlensing single-lens event with a 
measured Einstein radius Og < 10 yas, they noted that 
all four events had giant-star sources, which led them to 
initiate a long-term study to find all finite-source/point- 
lens (FSPL) events with such giant-star sources in 
the 2016-2019 Korea Microlensing Telescope Network 
(KMTNet, [Kim et_al]]2016) database. They were moti- 
vated in part by the fact that this was the second giant- 
source FSPL event from 2019 with 0g < 10 yas, the first 
being OGLE-2019-BLG-0551 (Mróz ot all 2020). They 
sought to place these FSPL FFPs in the context of a 
homogeneously selected sample that would include stars 
and brown dwarfs (BDs), as well as FFPs. For technical 
reasons that are briefly reviewed below, they were only 
able to probe the 2019 database at that time. Their sys- 
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tematic search yielded a total of 13 giant-source FSPL 
events, including the two FFPs. They carried out a 
variety of tests and concluded that there was nothing 
suspicious about the sample. However, they refrained 
from, and strongly cautioned against, drawing any sta- 
tistical conclusions from this sample, despite the fact 
that it was homogeneously selected. They argued that 
because their study was motivated by noticing two FFP 
detections in a single year, it could suffer from publica- 
tion bias. 

Two factors prevented (Kim et al], from carry- 
ing out a 4-year search immediately. First, the form 
of the search had been made possible by a recent up- 
grade to the 2019 online database, but this had not 
yet been extended to the previous 3 years. Second, for 
2019, had supplemented their auto- 
mated analysis of the online database with a special, 
more aggressive, execution of the EventFinder algorithm 
Kima al polsa) that was tailored to giant sources. A 
more aggressive search was required because the stan- 
dard EventFinder algorithm takes advantage of the fact 
that most single-lens/single-source (1L1S) events (and 
even most binary events) can be qualitatively matched 
to one of two variants of 2-parameter pro- 
files. However, this is not the case for typical FFP can- 
didates, for which the source radius, 0x, is often 20g, 
leading to a “top hat” light curve shape. These upgrades 
and new searches required more than a year of effort. 

The search proposed by [Kim et al] (2021) was there- 
fore carried out as these steps were completed season-by- 
season. In the course of this work, re- 
ported the discovery of KMT-2017-BLG-2820, the third 
FFP in the KMTNet sample, which by this time was 
the sixth FSPL FFP overall. That is, in the intervening 
time, had discovered another FFP, 
OGLE-2016-BLG-1928. pointed out 
that all six FFPs lay below an apparent “gap” in the 
cumulative 0g distribution from Figure 9 of [Kim ct al] 
(2021), which they dubbed the “Einstein Desert”. 

They argued from the definition of 65, 


On = 4/ KM Trel 


that such a desert would be a natural consequence of 
a bimodal function of the lens mass, M. That is, the 
main stellar and BD lens population in the bulge (with 
small lens-source relative parallaxes, mre) would con- 
tribute relatively small Einstein radii, but would be cut 
off at 0g ~ 30 yas due to a sharp drop in the BD mass 
function, e.g., near M ~ 0.02 Mo, combined with the 
paucity of phase space for 715 was (corresponding 
lens-source relative distances Dzs = Ds — Di <300 pc). 
Then, a second population of low-mass FFPs, e.g., 
M ~ 5 Mọ, with a much higher specific frequency than 
stars and BDs would give rise to the small-0p events. 
They then argued that a population of FFPs cannot 
be in the bulge and produce the observed distribution, 
because an analogous population in the disk would pro- 
duce larger-0g events that would “fill in” the gap. By 
contrast, if the small-Og events were explained by disk 


(1) 


FFPs, an analogous population in the bulge would not 
be detected. Specifically, for typical giant sources, the 
very numerous FF Ps in the bulge would not give rise to 
recognizable events because they would induce magnifi- 
cation, A, changes of only 


2 2KM Trel M Trel 0. 

£ = a 0.13 | — 

p? 62 5 Ma / \ 20 was / \ 6 pas 
(2) 


where p = 0,/0_21. Hence, the FFP events would 
mainly come from the less numerous disk population. 
Then, a high inferred specific frequency would be nec- 
essary to compensate for the overall relative paucity of 
disk lenses. (In Section [8] we will show that this reason- 
ing is only partly correct.) 


recognized that with only 13 FSPL 
events, the evidence from {Kim et al] for a desert 
was relatively weak. However, they adduced two addi- 
tional arguments. First, they noted that the apparent 
gap was not (yet) contradicted by the ongoing search for 
FSPL events in 2018 and (part of) 2017 data. Second, 
they pointed to an analogous desert in the distribution 
of Einstein timescales, 


A-1~ 


0 
tr = pE 
Hrel 


found by in their analysis of point- 
source/point-lens (PSPL) events extracted from 5 years 
of Optical Gravitational Lensing Experiment (OGLE) 
data. Here, 4e] is the lens-source relative proper mo- 
tion. The same argument predicts a desert in the dis- 
tribution of Einstein timescales, although it is expected 
to be less pronounced than the one in Einstein radii be- 
cause the former is a convolution of the proper-motion 
distribution with the latter. 


Thus, this picture gave a coherent account of all avail- 
able data at that time. Nevertheless, 
counseled a wait-and-see approach. If their hypothe- 
sis were correct, they predicted that the Einstein desert 
would remain relatively parched as new FSPL events 
were added. If not, one should expect that it would be 
gradually “filled in”. 


Mrel = [Hral ’ (3) 


Here, we present the full sample of 30 FSPL giant- 
source events from the KMTNet 2016-2019 database. 
We find that the Einstein Desert is substantially more 
distinct than in the 13-event sample from 2019. This 
strengthens the case for a bimodal mass function, with 
a broad peak of stars and BDs complemented by a large 
population of FFPs at much lower masses. 


We show that the FFPs are more numerous than the 
known bound planets. We raise the possibility that the 
heavier FFPs have mainly been ejected to very wide 
(rather than unbound) orbits, and may largely be ac- 
counted for by the still poorly traced population of very- 
wide-orbit planets. We also suggest that interstellar as- 
teroids may be the extreme end of the FFP population 
and, in particular, may follow the same power law. 


S 
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2. REVIEW OF DETECTION PROCEDURES 


The FSPL events reported here are derived from the 
KMTNet microlensing survey, which uses three identi- 
cal 1.6m telescopes, each equipped with 2° x 2° cam- 
eras (2016). The telescopes are located in 
Chile (KMTC), South Africa (KMTS), and Australia 
(KMTA). See Figure 12 from|Kim et al. for the 
field locations and cadences. 

We essentially follow the same procedures for select- 


ing FSPL preliminary candidates as was described by 
(2021). We refer the reader to that paper 
for details. Here we give only a brief overview, and we 
also describe two important differences between their 
2019 sample and the 2016-2018 additional events that 
we present here. One of these was already mentioned 
by (2021), while the other is new. 

We first conduct a purely automated search of the on- 
line summary table of KMT events (about 3000 events 
per year), making use of the impact parameter uo, the 
Einstein timescale tg, and the source magnitude Ig 
(from the pipeline Bacarisklosg psu fit), together 
with the tabulated field extinctior}] Az, to impose two 
selection conditions, 


Is.0,est =TIg— Ay < 16, (4) 
and 
0. est =ï 
Hthresh = — > I masyr `; 
UotE 
bz est = 3 pas x 1006708 -41))/5, (5) 


The first condition restricts the sample to giants. On 
rare occasions, it allows foreground main-sequence or 
subgiant stars, but these are eliminated at the stage of 
detailed investigations. 

The second condition rejects the great majority of 
PSPL events while accepting (for further investigation) 
the overwhelming majority of FSPL events that satisfy 
Equation (4). First note that the 0, est estimate implic- 
itly assumes that the source has a similar color as the 
clump. This will be approximately so except for stars 
well up on the giant branch, for which it will be an 
underestimate. Then, assuming that the online fits for 
PSPL events (i.e., p < uo) are approximately accurate, 
PSPL events will have 


-1 -1 
ype?! — 0.33 ZES (vo tg Os est 
thresh. =" yr (1/3 20 day 6 uas j ` 
Hence, most PSPL giant-source events will fail this crite- 
rion, leaving a tractable subset to be rejected by manual 
light-curve fitting. 

On the other hand, for FSPL events, we do not expect 
the automated fit to yield the correct uo. For p < 1, 
we expect UO fittE, fit = teff fit = tx true = PtruelE,true to 
the extent that the fit is influenced by the width of the 
peak, and 1/uo,fit > Amax, fit > Amax,true > 2/Ptrue to 
the extent that the fit is influenced by the height of the 


+The KMTNet webpage adopts: Ay = 7 Ag, where Ax is from 
2012. 


peak. These imply Hrel ~ Hthresh and Hre ~ Lthresh /2, 
respectively. Hence, FSPL events that fail Equation (5) 
have relative proper motions [rel masyr—!. Using 
the approach in the Appendix to [ould at all (2021), 
the probability for a lens to have such a small proper 
motion is p < (tUthresh/Ou)?/6V7 —> 4 x 107%, where 
we have adopted an isotropic bulge proper-motion dis- 
persion of op ~ 2.9masyr~!. Table 4 of [Kim et al] 
(2021), which lists both the estimator wthresh and the 
true value Hre] for 13 events shows that the former is 
conservative in the sense that plrei/{Uthresh Sl. Only two 
events (OGLE-2019-BLG-0551 and OGLE-2019-BLG- 
1143) violate this inequality, and they do so only mildly, 
both with ratios 1.2. Note that neither of the two 
FFPs in this sample (OGLE-2019-BLG-0551, p = 4.5; 
KMT-2019-BLG-2073, p = 1.1) are in the p < 1 
regime, but both approximately satisfy the inequality, 
with Hrel/Hthresh = 1.2 and 0.6, respectively. 

Each such candidate, typically of order 60 per year, is 
first reviewed visually. Of order 10% are eliminated at 
this stage for various reasons, primarily that they are not 
single lens events or that there are no data near the peak, 
which would be required to measure p. The remainder 
are fitted to PSPL and FSPL forms using pipeline KMT 
data. If FSPL is preferred by Axĝpı > 3 and the best 
fit has a normalized impact parameter zo = uo/p < 1 
then the data are subjected to tender loving care (TLC) 
re-reduction. Then, if AXĝspl > 20, the event is selected 
as FSPL. For a few cases, the pipeline light curve has 
an obvious FSPL form for which PSPL would give a 
very poor fit. For these, we move directly from inspec- 
tion to TLC. Note that both the pipeline reductions and 
TLC reductions use pySIS {Albrow et all Roog), which 
is a specific implementation of difference image analysis 

For 2019, separately searched 
the output of the real-time event selection algorithm 
AlertFinder (Kim et all, and the end-of-year se- 
lection algorithm EventFinder (Kim ct all, 2018A). Be- 
cause AlertFinder did not operate in the wings of the 
season and also because it searches only for rising events, 
there were wing-of-season events that it missed but were 
found by EventFinder. However, even in the region of 
overlap, there were FSPL events that were missed by 
one or the other. For example the FFP FSPL event 
OGLE-2019-BLG-0551 was found by AlertFinder but 
not EventFinder. These discrepancies motivated the ad- 
dition of an aggressive giant-source search, which had 
two purposes. First, it served as a check that the union 
of the AlertFinder and EventFinder samples for 2019 
was effectively complete. And, indeed, the special giant- 
source search did not yield any additional FSPL events. 
However, [Kim et al] also noted that 2019 was the 
first year for which AlertFinder was in full operation. 
For 2018, AlertFinder operated only beginning in June 
and only in the northern bulge. It did not operate at all 
in 2016 and 2017. Hence, given that EventFinder was 
missing FSPL events in 2019, the special giant-source 
search would have to be applied to previous seasons. See 


(2021)) for further refinements of the giant- 
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source search. 

Finally, we modified the search procedures relative 
to better deal with saturated data. In 
2019, events that were saturated at peak were simply 
eliminated as uninterpretable at the visual-inspection 
stage. However, we subsequently realized that many 
such events could be analyzed by incorporating the V- 
band data for the regions of the light curve that were 
saturated in the J band. For example, given the typical 
intrinsic colors of giants (V — I)o ~ 1.1, typical redden- 
ing of KMT fields E(V — I) ~ 1.8 and the roughly 0.65 
deeper photometric zero point in the V band, the on- 
set of saturation would be ~ (1.1 + 1.8 — 0.65) = 2.25 
mag fainter in the V band than the J band. Of course, 
V-band data were taken 10 times less frequently in 2017- 
2019 than J band, and even less frequently in 2016. 
However, for events that are well characterized by bright 
(so, generally small-error) -band light curves, only a 
few points are needed over peak to measure p. We ap- 
plied this approach to the saturated events from 2019, 
but this did not lead to new FSPL events. However, 
it did prove useful in fully characterizing FSPL events 
from 2016-2018. 


3. EXPECTED g DISTRIBUTION 


Our primary motivation is to probe the poorly under- 
stood population of FFPs. To do so, it is useful to 
understand the 0g distribution that is expected from 
stars and BDs. We first show that a substantial ma- 
jority of these detections will come from bulge lenses 
and second that these are expected to cover the range 
30 pasSOp S300 pas. We then turn to the role of disk 
lenses. 

Ignoring selection effects, the rate of microlensing 
events for a population of fixed mass M and toward 
a given source is 


Tan x J tscthe( D1) DEAD ivf iter)tuerg( Ms Dr) (7) 


where v = v( Dz) is the density at distance Dz, f(u) = 
f (frei; DL) is the 1-dimensional proper motion distri- 
bution and g(M, Dr) is the mass function. However, 
for FSPL events, we replace 0g(Dr) > pOn(Dr) = 0x, 
which eliminates one instance of implicit dependence on 
Dy, (and M; [Gould & Yee)2013). We note that with the 
exception of very nearby lenses D;,<1 kpc (see below), 
the mean proper motion (Hre1) is approximately inde- 
pendent of distance, while we also approximate g(M) as 
being independent of distance. Hence, for FSPL events, 


Ds 
T FSPL X 8. (soa) 9(M) | dD, D?v(Dr). (8) 
D,=0 


That is, [pgpz is proportional to the total column of 
lenses in the observation cone between the observer and 
the bulge source. Without any detailed model, it is clear 
that this column is dominated by bulge lenses: other- 
wise, e.g., there would be a cloud of disk clump stars 
trailing toward brighter magnitudes from the observed 


bulge clump on the color-magnitude diagram (CMD). 
Hence, bulge stars will dominate most of the 0g distri- 
bution, while disk lenses will dominate only in regions 
that are inaccessible to bulge lenses. 

The bulge mass function is well populated over the 
range 0.02 MoSM<1 Mo. The upper limit is set by 
the fact that the bulge is primarily an old popula- 
tion for which the stars that were born M21 Mo 
have now mainly evolved to become remnants, mostly 
white dwarfs. The lower limit is an approximation for 
the steeply falling mass function in the BD regime. 
The bulge has a depth of order 1 kpc, implying that 
the typical range of lens-source relative parallaxes is 
0.005 masS7¢1 {£0.02 mas, corresponding to lens-source 
separations 0.35 kpc<(Ds — Dr)<1.5kpc. Above the 
upper limit, the product of the disk and bulge density 
distribution declines sharply. Below the lower limit, 
the amount of available phase space declines. None 
of the four limits just described is a hard boundary, 
but together they indicate a well-populated range, for 
Op = (KM Trel of 30 pasSOp<400 pas. 

Then, we should consider whether these intrinsic 
boundaries are impacted by additional detection effects. 
For typical giant-star source radii 0, ~ 6 uas, these lim- 
its correspond to Amax ~ 2/p ~ 10 and 135, respec- 
tively. The lower limit is not near any selection bound- 
ary, but saturation becomes an issue at the upper limit. 
For typical giant sources Jy ~ 14.5 and extinctions A; ~ 
2, the upper boundary implies Ipeak ~ 11.2, which is well 
into the range of KMT saturation, whose onset is seeing 
(and hence observatory) dependent. As we discussed in 
Section [2] it is sometimes possible to recover from sat- 
uration using V-band data. Nevertheless, the onset of 
saturation steepens the decline in the upper range of 0g 
detections, leading us to adopt 30 pas{<0g 300 pas as 
the expected range for bulge lenses. 

If we were to adopt a similar mass function for disk 
lenses, then they would populate a range that is trans- 
lated upward by a factor \/7rel,disk/Trel,bulge ~ 2 in Op, 
where we have adopted Tyel,disk ~ 60 was and Trel,bulge ~ 
15 pas as “typical” values for the disk and bulge, respec- 
tively. Hence, the disk lenses would appear as “sprinkled 
in” (and not generally individually identifiable) over the 
range 60 wasS6¢<400 uas, and they would only start to 
dominate for 6g2400 pas. The fact that the disk mass 
function extends above MZ Mo augments this effect but 
does not qualitatively change it. As just mentioned, this 
high 0g regime is suppressed by saturation effects. How- 
ever, another way to express the relation of the disk and 
bulge contributions is that at fixed 0g, the disk lenses 
have ~ 4 times lower mass than the bulge lenses. Hence, 
at 0g ~ 400 yas, they are near the peak of the mass func- 
tion M ~ 0.3 Mo, whereas the bulge lens mass function 
at M ~ 1 Mo is much lower. Hence, we expect that, 
contrary to the main part of the 0g distribution, the 
upper range will be a comparable mix of disk and bulge 
lenses. 

We note that the higher proper motions of nearby 
lenses Dz{1kpec do not compensate for their low 
observing-cone volume in the rate equation. The mean 


Name 
OB180705 


KB180244 


OB180626 


MB17147 


OB171254 


OB170560 


OB170905 


MB17241 


OB170084 


OB161045 


KB161128 


OB161540 


KB162057 


MB16258 
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MICROLENS PARAMETERS FOR FSPL GIANT-STAR EVENTS 


KMT Name 
KB181882 
(errors) 
KB180244 
(errors) 
KB182309 
(errors) 
KB170132 
(errors) 
KB170374 
(errors) 
KB172830 
(errors) 
KB171022 
(errors) 
KB170818 
(errors) 
KB170726 
(errors) 
KB160848 
(errors) 
KB161128 
(errors) 
KB162262 
(errors) 
KB162057 
(errors) 
KB160606 
(errors) 


Table 1 
to uo 
8316.30592 0.00123 
0.00145 0.01239 
8312.42231 0.27149 
0.00118 0.00301 
8230.28800 0.16989 
0.00132 0.00255 
7850.99388 0.09138 
0.00043 0.00097 
7952.25224 0.00000 
0.00025 0.00000 
7859.52120 0.02900 
0.00100 0.01800 
7895.73816 0.08702 
0.00401 0.00221 
7883.47648 0.21697 
0.00057 0.00437 
7807.13435 0.01226 
0.00212 0.00022 
7559.20148 0.01061 
0.00115 0.00208 
7486.32769 0.00951 
0.00147 0.00105 
7606.72400 0.60500 
0.00100 0.02700 
7467.93171 0.00000 
0.00196 0.00000 
7537.29800 0.18800 
0.00100 0.00400 


tr 
46.247 
0.147 
4.430 
0.027 
2.983 
0.024 
2.695 
0.016 
15.278 
0.046 
0.901 
0.002 
7.734 
0.046 
1.845 
0.024 
43.642 
0.501 
12.030 
0.099 
12.550 
0.428 
0.330 
0.003 
11.374 
0.105 
3.722 
0.045 


P 
0.01151 
0.00503 
0.35487 
0.00283 
0.21484 
0.00251 
0.13451 
0.00102 
0.02527 
0.00009 
0.89000 
0.00300 
0.15031 
0.00254 
0.29496 
0.00483 
0.02335 
0.00028 
0.02943 
0.00191 
0.02586 
0.00105 
1.63100 
0.00800 
0.06645 
0.00078 
0.57400 
0.00300 


fs,KMTC 
5.5354 
0.0187 
11.0566 
0.1300 
10.3031 
0.1520 
4.3062 
0.0393 
0.7068 
0.0024 
18.5140 
0.0010 
83.4744 
0.9971 
4.8304 
0.1019 
0.4049 
0.0055 
1.3546 
0.0141 
0.2703 
0.0106 
19.6960 
0.0000 
2.4595 
0.0309 
49.7270 
0.4670 


Notes: The units of to and tg are HJD’ = HJD — 2450000 and days, respectively. Fluxes are in units of an 
names are abbreviations for, e.g., OGLE-2018-BLG-0705. For KB172820, selaa] 
Shvartzvald et al] (2019). For OB171186 (=KB170357), see|Li etal] (2019). For OB180705, the fit included parallax: (Tmp, N, Tp, E) = 


(—0.324, —0.017) + (0.139, 0.008). 


Name KMT Name 


OB180705 KB181882 
KB180244 KB180244 
OB180626 KB182309 
MB17147 ~=KB170132 
KB172820 KB172820 
OB171254 KB170374 


OB170896 KB170799 
OB170560 KB172830 
OB170905 KB171022 
MB17241 KB170818 
OB171186 KB170357 
OB170084 KB170726 
OB161045 KB160848 
KB161128 KB161128 
OB161540 KB162262 
KB162057 KB162057 


On 
1285.00 
96.72 
29.00 
46.24 
5.94 
224.11 
139.60 
34.90 
256.15 
25.93 
93.90 
204.11 
249.31 
212.79 
8.80 
109.71 


I = 18 system. Event 


(2021). For OB170896 (=KB170799), see 


Table 2 
CMD AND DERIVED PARAMETERS FOR FSPL GIANT-STAR EVENTS 
(V — I)o Io Os 
1.46 13.14 14.78 
2.26 12.06 34.25 
1.20 14.69 6.53 
0.94 14.18 6.22 
1.13 14.31 7.05 
1.27 14.97 5.68 
1.26 14.95 5.71 
2.43 12.47 31.10 
1.60 11.19 38.50 
1.24 14,29 7.65 
2.12 12.39 26.90 
1.33 15.41 4.78 
1.23 14.37 7.33 
1.00 14.43 5.49 
1.70 13.58 14.36 
1.20 14.58 7.29 
2.92 12.31 43.89 


MB16258 KB160606 


Notes: The units of 6, and 0g are yas, while those of Hre} and Lthresh are mas yr™-. 


KMT-2018-BLG-1882, OGLE-2018-BLG-0705, and MOA-2017-BLG-147 


76.45 


Hrel Hthresh 20 = uo/p 
10.10 13.31 0.104 
7.97 7.79 0.764 
3.55 5.63 0.791 
6.26 13.75 0.679 
7.95 10.34 0.150 
5.35 10.14 0.000 
3.42 5.94 0.115 
14.17 5.80 0.033 
12.09 5.49 0.579 
5.13 5.09 0.736 
2.65 2.88 0.319 
1.71 2.38 0.525 
7.55 43.72 0.360 
6.19 24.10 0.368 
9.74 8.80 0.371 
3.52 8.41 0.000 
7.50 1.65 0.328 


1 


Event names are abbreviations for, e.g., 
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proper-motion term in this regime scales (Hreal(Dr)) ~ 
V1 typ/ Dr where v1 typ ~ 40kms~! ~ 8kpcmasyr~? is 
the typical transverse velocity of nearby stars. Hence, 
when the proper-motion term is moved inside the in- 
tegral, it takes the form v1 typ f dDr Dz... in place 
of (mra) f dDr D? ... at larger distances. That is, the 
nearby lens contribution is still heavily suppressed by 
the volume factor. 

Finally, it is reasonable to expect sensitivity to FSPL 
events with peak magnifications down to Amax ~ 2, i.e., 
p ~ 1.4 (Equation (2), hence 0g ~ 0x/1.4 ~ 4 uas. Thus, 
the survey should have sensitivity that goes almost a 
decade below the Einstein radii that are generated by 
known stars and BDs. 


4. 30 FSPL EVENTS FROM 2016-2019 


Table [shows the microlens fit parameters for the 14 of 
the 17 new FSPL events that were discovered in 2016- 
2018 data. The parameters for the remaining three are 
adopted from previous publications. See Appendix [A] 
For the 13 FSPL events from 2019, see Table 3 of 
(2021). 


For the color-magnitude analysis, we used the same 
procedures (2021), ultimately based on 
(2004). Briefly, we perform a special reduc- 
tion using pyDIA i017), which carries out 
light-curve photometry and field-star photometry on the 
same system. We find Is by regression of the J-band 
light curve on the best model and the (V — I)s color 
by regression of the V-band on the J-band light curve. 
We find the clump centroid of the field stars, and so 
calculate the offset A[(V — I), I] of the source relative 
to the clump. We estimate the intrinsic clump color 
(V — Dao = 1.06 from (2013) and its 
intrinsic magnitude «ao from Table 1 of 
(2013). We convert from oe to V/K using the color- 
color relations of (1988). Finally, we 
estimate 0, using the color/surface-brightness relations 
of [Kervella et al! for K giants and [Groenewegen] 
(2004) for M giants. The results of the CMD analysis are 
shown in Table P] for the 17 events from 2016-2018. For 
the 13 FSPL events from 2019, see Table 4 of[Kim ct al] 


(2021). 


Among the 17 newly reported FSPL events, we re- 
covered the FFP OGLE 2016 BLG-1540 iM etal 
2019) (and, as already reported by 
the FFP KMT-2017-BLG-2820), but did not find any 
additional FFPs. Hence, there are a total of four 
FFPs in the sample of 30 FSPL events, OGLE-2016- 
BLG-1540, KMT-2017-BLG-2820, OGLE-2019-BLG- 
0551, and KMT-2019-BLG-2073. 

We did not recover either of the other two known 
FSPL FFPs, OGLE-2012-BLG-1323 (Mróz ct al), [2019) 
and OGLE-2016-BLG-1928 (Mróz et all, |2020). The 
first preceded the KMT survey. For the second, the 
event was not found either by EventFinder nor in the 
special giant-source search because there are only 4 mag- 
nified KMT points. Even if identified, it could not 
have been reliably characterized from these four points. 
OGLE?’s original identification of this event was post- 


season and was based on 10 magnified points that probe 
both the peak and the wings. The role of KMT data (in 
addition to confirming the event from KMTC) was to 
rule out binary models by the flat behavior of KMTS 
data starting about 8 hours after the last OGLE point. 

In Appendix[A] we provide notes on individual events. 


5. SAMPLE CHARACTERISTICS 


Notwithstanding the justified caution of 
(20211), the results from 2019 prove to be reasonably rep- 
resentative of the four year sample. For the four years 
2016-2019, the FFP/FSPL counts are (1/5, 1/9, 0/3, 
2/13). For 2016-2018, they are 2/17 compared to 2/13 
for 2019, which is consistent with normal Poisson varia- 
tions. If four detections are distributed randomly among 
four seasons, the most common configuration (144/256) 
is the observed one, (2,1,1,0). The one feature that is 
possibly unusual is that there were only 3 FSPL events 
in 2018 compared to 13 in 2019, so one may wonder 
whether the seasonal distribution of FSPL events, n; = 
(5,9, 3, 13) can also be considered as a “typical” outcome 
of Poisson sampling. We investigate this by adopting an 
annual expectatior<] A = 30/4 = 7.5 and construct a log 
likelihood InL = X$; In(e~*A"* /n;!) = —10.31. We 
can then compare this to the same statistic for 10° re- 
alizations in which 30 events are randomly distributed 
among 4 seasons. We find that 7.8% have smaller like- 
lihood. Thus, this outcome is not particularly “unex- 
pected”, especially given that this is a posterior test. 
Therefore, neither the seasonal distribution of the FFPs 
nor that of the FSPL events can be considered unusual. 
presented several diagnostic plots 
whose purpose was to probe for statistical artifacts 
and/or patterns in the sample and to identify individual 
outliers that might require explanation and/or further 
investigation. We present analogous plots here. 

Figure [I] shows the color-magnitude diagram (CMD) 
source positions of the FSPL events relative to the 
clump, color-coded by season. The blue (2019) points 
and the black background points (taken from the dered- 
dened CMD of KMT-2019-BLG-2555) are the same as 
those shown in Figure 8 of {Kim et al] (2021). Of the 30 
sources, 15 are clump stars (or giant-branch stars that 
are superposed on the clump), 4 are lower-giant-branch 
stars, 9 track the upper giant branch quite closely, one 
(KMT-2019-BLG-1143, blue) lies significantly below the 
upper giant branch, and one (OGLE-2017-BLG-0905, 
green) lies well above the upper giant branch. 

The broad pattern of the FSPL events relative to 
the background is as expected: they are concentrated 
in the clump (i.e., the most densely populated region 
of giant stars), and are somewhat over-represented on 
the upper giant branch and under-represented on the 
lower giant branch, which is expected due to larger and 


2Tnspection of the formula for In L shows that the choice of À plays 
no role: it just adds a constant 30 1n A— 4A — 30.45 to each trial. 
That is, we could have just as well used ln L = — Dai Inni!. 
We keep the more complex form to maintain familiarity for the 
reader. 
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Figure 1. Color-magnitude locations of the 30 FSPL events 
relative to the clump, shown against a background of field 
stars from the CMD of KMT-2019-BLG-2555, i.e., the same 
as was used in Figure 8 of fozi). The events are 
color-coded by discovery year. FFPs have red cores. As ex- 
pected, the sources are concentrated in the clump (the high- 
est concentration of giants), are somewhat over-represented 
on the upper giant branch (due to larger cross sections), and, 
conversely somewhat under-represented on the lower giant 
branch. The largest outlier, OGLE-2017-BLG-0905 (green 
point near top), is discussed in Appendix [A] while the sec- 
ond largest, KMT-2019-BLG-1143 (blue point toward the 
right), was discussed by [Kim et al] aa The lack of de- 


tections with ordinate >15.5 is discussed in Section [5] with 


reference to[Kim et al] (2021). 


Figure 2. Field positions of 30 FSPL events in Galactic coor- 
dinates, color-coded accorded to year of discovery (see labels 
in Figure {I}. The background black points are (as in Fig- 
ure 4 of[Kim et _al][2021) EventFinder events from 2019. As 
in that Figure, the KMT fields are shown as black squares, 
with blue field numbers. Their observational cadences can 


be obtained from Figure 12 of[Kim et al] (2018al). 


30 


20 
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Figure 3. Cumulative distribution of normalized impact pa- 
rameters zo = uo/p. The red line connects the first and last 
points, while the green line connects the first and penulti- 
mate points. The underlying distribution is rigorously uni- 
form in zo, so the eye is struck by the absence of detections 
toward zo œ~ 1. According to a Kolmogoroyv-Smirnov (KS) 
test, this is not significant, but the effect may be real because 
the KS test does not capture all the physics. See Section [5] 
The distribution relative to the green line suggests that the 
sample is homogeneously selected for zo<0.8. 


smaller cross sections, respectively. The complete ab- 
sence of FSPL events between 16 < I < 15.5 would 
appear to be somewhat surprising. However, {Kim et al] 
showed that this is due to systematic underes- 
timation (with considerable scatter) of source flux for 
FSPL events by the automated pipeline. See their Fig- 
ure 7. 

The outlier KMT-2019-BLG-1143 was investigated 
by [Kim et al] (2021). We investigate the other outlier 
(OGLE-2017-BLG-0905) in Appendix [A] 

Figure] shows the distribution of the 30 FSPL events 
on the sky in Galactic coordinates, against the same 
background of 2019 EventFinder events that appears in 
Figure 4 of [Kim et al] (2021). The FSPL events are 
concentrated toward the Galactic plane, as would be 
expected for giant-source events because they are more 
heavily represented in regions of high extinction. How- 
ever, this is a weak effect, and no strong statistical con- 
clusion can be drawn from this apparent pattern. 

Figure [3] shows the cumulative distribution of the 30 
FSPL events as a function of normalized impact param- 
eter zo = uo/p. In nature, this distribution is exactly 
uniform, so any (statistically significant) deviation from 
a uniform distribution should reflect selection effects. 
had argued based on their analogous 
Figure 3, that the distribution was consistent with being 
uniform over [0 < zp < 1]. Comparison of the observed 
distribution with a uniform expectation (red line) in- 
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Figure 4. Cumulative Distribution of log@m for 30 FSPL 
events with giant-star sources that were found in four sea- 
sons of KMTNet data, 2016-2019. The absence of detections 
On<4 pas is due to selection effects, while the paucity for 
022300 pas is due to a combination of selection effects and 
the lack of bulge stars M21Mo. The remaining feature, 
the factor ~ 3 “Einstein Desert”, 8.8 pas < 0m < 26 pas, 
is the signature of two populations, i.e., BDs/stars with 
masses 0.02<M/Mo<1 to the right and FFPs with much 
lower masses to the left. The former are detected mainly 
in the bulge, while the latter are likely to be about equally 


drawn from the bulge and the disk. See Section [8.3] This 
Figure can be compared to Figure 10 from|Kim et al} (2021). 


dicates that, with accumulating statistics, this may no 
longer be true. A Kolmogorov-Smirnov (KS) test yields 
an unimpressive p = 28%. However, KS is a very “for- 
giving” test because it assumes no prior information on 
the form of the deviation. By contrast, our eye is drawn 
to the near absence of events with zo > 0.82, where we 
might expect the greatest difficulty for detections. The 
green line, which ignores the single detection at zo œ 1, 
appears more satisfying. argued that 
there is substantial information about p in the regions 
of the light curve having z = u/p > 1, so that uniform 
sensitivity over [0 < zo < 1] was plausible. However, it 
is also plausible that there would be at least some effect 
for z ~ 1, as seems to be indicated by Figure [B] albeit 
weakly. 


One might suspect that the dearth of events with 
zo ~ 0.9 is the result of the z9 < 1 boundary in our ini- 
tial selection combined with ordinary statistical noise. 
That is, some events with initial estimates zp init ~ 0.9 
were ultimately eliminated after they were determined 
to have zo,rtc > 1 from subsequent TLC reductions, 
but the “complementary” events that would have had 
zo,rLc ~ 0.9 (if TLC reductions had been done) were 
not investigated after it was found that zo init > 1. In 
fact, we were concerned about this and obtained TLC 
for 5 events with zo jinit21. However, zo,rLc > 1 was 


confirmed for all 5. Hence, we do not believe that such 
selection bias is a strong effect. 


6. THE EINSTEIN DESERT 


Figure[4]shows the cumulative distribution log ôg for the 
30 FSPL events in the survey. It exhibits three principal 
features; 1: a paucity of detections with 6_2300 „as, 2: 
a complete absence of detections with 0g < 4 as, and 
3: the “Einstein Desert”, i.e., the absence of detections 
over the interval, 8.8 was < 0g < 26 Op. 

As discussed in Section B] the first was expected, pri- 
marily due to the rapid fall-off of the bulge mass function 
for M21 Mo, but somewhat augmented by selection ef- 
fects due to saturation. The latter will be examined 
more closely in Section [7] 

As will also be discussed in Section [7] the second fea- 
ture is likewise due to selection. 

However, selection effects play no role in the third fea- 
ture, the “Einstein Desert”, because the selection func- 
tion is smoothly increasing from a factor of 2 below the 
Desert to a factor of a few above it. 

As was already anticipated in Section [B] the upper 
shore of this desert (at 0g = 26 yas) can be attributed 
to the sharp fall-off in the bulge mass function in the BD 
regime. In this sense, it is qualitatively similar to the 
first feature, but with a crucial difference. Being gener- 
ated by a fall-off at the high end of the mass function, 
the first feature is “washed out” by contributions of disk 
lenses of much lower mass (hence, higher specific fre- 
quency) but with similar 0g = VKM Tre- By contrast, 
the population of disk stars and BDs only contribute to 
Og in regions well above the desert, and so the upper 
edge of the Desert is sharp. 

The 4 events that lie below the Desert must repre- 
sent a separate, low-mass population because they have 
Einstein radii that are a factor 23 below those of the 
lowest-mass bulge BD lenses (and even farther below 
those of the lowest-mass disk BD lenses). The existence 
of the Desert is a powerful constraint on the nature of 
this population: whatever model one adopts must not 
only reproduce the observed low-0g detections but must 
also respect the absence of such detections in the Desert. 

Figure]shows a scatter plot of FSPL proper motions 
against Einstein radii. Generally, the proper motions are 
consistent with those expected for bulge and disk lenses. 
Note that the FFPs are within this range, somewhat 
more tightly grouped but with a median similar to the 
sample as a whole. 


7. SELECTION EFFECTS 
7.1. Selection Effects Due to Lens Mass (thus, 0g) 


There are two principal selection effects due to lens mass 
that can prevent a given event from entering the sam- 
ple. These derive from the requirements that, first, the 
light curve must initially be selected as a “microlens- 
ing event”, and second, that finite source effects must 
be detected in this event. In Section 2] we adopted a 
threshold for finite source effects: AxXf&pı = X? (PSPL) — 


y?(FSPL) > 20. The threshold for event detection, 
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Figure 5. Scatter plot of log Hre, vs. log 0m for the 30 FSPL 
giant-source events from four years of KMTNet data. The 
events are color coded by year: cyan (2016), green (2017), 
magenta (2018), and blue (2019). The four FFPs at the 
left (with red interiors) all have pire) > 4mas yr’, whereas 
the 26 non-FFPs are distributed more broadly. However no 
strong conclusion can be drawn from this. Note that Figure 9 
of (20211) corresponds to the blue points in this 


Figure. 
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Figure 6. Measurements (circles) and sensitivity limits (error 
bars) of 0g for 29 FSPL events, rank ordered by lower limits. 
As the 0g of hypothetical events is progressively reduced 
below that of the actual event, it eventually falls below the 
detection threshold, either because Ay2.1.¢, < 1000 (required 
to detect the event) or AX%pı < 20 (required to be chosen 
as an FSPL event). As ĝm is progressively increased, the 
peak of the event eventually saturates, making it impossible 
to measure Oz. KMT-2019-BLG-2528 with Ax2.1.c, = 973 is 
excluded because the detection itself fails AXZ lect > 1000. 
The magenta histogram shows the number of FSPL events 
that were sensitive to each On. 


AX? jects Łe, the difference in x? between microlensing 
and non-microlensing interpretations of the light curve 
as determined by EventFinder and the special giant-star 
searches, varied somewhat for different subsamples (see 
[Kim et al][2021). For purposes of this section, we im- 
pose the minimum threshold that was common to all 
subsamples, namely, Ax2.),., > 1000. This has the ef 
fect of eliminating one of the FSPL events from the sam- 
ple, i.e., KMT-2019-BLG-2528 with AXŽject = 973. A 
hypothetical event can drop out of the sample either be- 
cause AXĝpı falls below our threshold of 20 or because 


Ax2 ject falls below our threshold of 1000. Under our 
assumptions for real versus hypothetical events, we will 
see that both AXĝ&pı and Ax2.1¢¢, Will monotonically de- 
cline with falling mass. 

To understand how selection effects impact these two 
criteria, consider two events that are “identical” in every 
respect except for the lens mass, e.g., a real event and 
a hypothetical event with either larger or smaller lens 
mass. Here, “identical” means that the lens and source 
of the hypothetical event both have the same 6 phase- 
space coordinates as the real event and that the source 
radius and temperature are also the same. Under these 
assumptions, flrel, Zo = Uo/p, tx = pte, and tef = uote 
are the same for the two events. However, because 0g 
scales as M'/?, tp, uo, and p all differ. Specifically, 
te x M/?, uo x M~'/?, and p x M~1/2. We also 
assume that the observational sequences of the real and 
hypothetical events are the same. Then, we can eval- 
uate upper and lower limits on @g for each real event 
based on comparison to hypothetical events under these 
assumptions. 

First consider the case that the real and hypotheti- 
cal events satisfy preal < Phyp < 1, i.e., in particular, 
Mhyp < Mreai. Then, because z is the same in both 
cases, both the PSPL and FSPL models will have essen- 
tially the same morphologies in the peak region of each 
event. In addition, because the duration of the tran- 
sit, Attrans = 2te y p? — uå, is invariant, they will have 
roughly the same number of data points over the peak 
region. On the other hand, they will have a magnifica- 
tion ratio Ahyp/ Areal > (Pnyp/Prea1) ~} = OE hyp/0E, real- 
Hence, the real event will be brighter than the hypothet- 
ical one, and the number of “signal photons” (difference 
between the true number for the FSPL model and the 
expected number for the incorrect PSPL model) will fall 
linearly with Of hyp, while the number of background 
photons remains the same. Hence, Ay?,,, will fall lin- 
early with Op nyp if the hypothetical event is “above sky” 
and quadratically if it is “below sky”. Once 0p falls suf- 
ficiently that p21, the functional forms of the models 
in the two cases start to differ, but the principle re- 
mains the same. We estimate Ay?,,, by scaling to the 
actual value found from the FSPL and PSPL fits to the 
real data and calculating the signal-to-noise ratio of each 
case by assuming typical 2” FWHM seeing and a typical 
I-band background of 19.25 mag per square arcsec. 

Next we turn to Ay?.,.,. As long as pS1, the prin- 
cipal impact of declining ôg on Ay?,,., is that tp is 
shorter, so that there are fewer magnified data points 


10 ANDREW GOULD 


for a given event. Hence, Ax2,),, declines as AXZ lect X 
te x Og. Eventually, 0m becomes small enough that 
p 2, 1, at which point, finite source effects dominate the 
morphology of the light curve. At this point, the effect 
of the further decline of 0m through p = 1 on AXŽyject is 
similar to that on AXĝpı: 


On the other hand, as the mass is increased, eventu- 
ally the peak flux freak ~ 2fs/p x M'/? will grow so 
large that the images become saturated over peak. As 
discussed above, in many cases it will still be possible 
to measure p from the V-band images. However, these 
are 10 times less frequent, and so for low-cadence fields 
they may not sample the peak region. Moreover, as the 
mass continues to be raised, eventually the V-band im- 
ages will also become saturated. We model this process 
by assuming the same seeing and background as above, 
and by reducing Ax? by a factor 20 for the regime in 
which J is saturated but V is not. Here, a factor 10 
comes from the lower cadence in V band and 2 comes 
from the lower photon counts in the V band. This sets 
an upper limit on the detectable 0g for a given event. 


In most cases, we do not actually know the lens mass, 
but we do know @, and 0g. Hence, we can still com- 
pare the real and hypothetical events in terms of a mea- 
surable quantity, 62 = 0,/p. Then, for each event i, 
with measured Einstein radius, 6p,;, there is some range 
OE,min,i < Ôe < OB max,i over which it could have been 
detected, where 0%, min,i < Ôg, is set by the Ay? thresh- 
old (either AXĝspı or AxZetect) and On,max,i > Ôp, is set 
by saturation. 


Figure [6] shows these ranges of sensitivity to 0m for 
the 29 FSPL events with Ay?2,),., > 1000, rank ordered 
by the 0% min of each event. The most important feature 
of this Figure is that the 4 lowest-6p events (i.e., the 4 
FFPs) are all pressed up near their detection limit. By 
contrast, and with one very telling exception, there is 
no such effect for the highest-@g events. For example, 
the 8 events with 0g > 200 uas occupy a broad range of 
positions relative to the detection limits for that event. 
This suggests that the fact that there is only one event 
with 0g > 500 pas is due to the low frequency of such 
events in nature rather than lack of sensitivity of the 
survey. Indeed, as shown by the magenta curve, of the 
29 FSPL events, 22 would be able to detect events with 
0g = 500 pas. Moreover, there is a well-understood rea- 
son for the paucity of 62500 uas events: as discussed 
in Section [B] these are almost entirely due to disk lenses, 
and in a mass range for which the mass function is al- 
ready in decline. Indeed, the one event that lies above 
this limit (OGLE-2018-BLG-0705), is known to lie far 
in the foreground Dz ~ 2.5 kpc by two independent ar- 
guments. First, it has a measured 7g (in addition to 
the known 6g), yielding mre, = Ogre. Second, its Gaia 
proper motion indicates that it is a member of NGC 
6544, whose distance is well-determined by several tech- 
niques. See Appendix[A] Hence, the fact that the low-05 
events are all pressed up against the detection limit hints 
that these may represent the “edge” of a population that 
we are just barely detecting. 


7.2. Selection Effects Due to rel 


There are also selection effects due to lens-source 
relative proper motion. To understand these, we 
consider hypothetical events with the same g and 
lens-source trajectory, but moving with faster or slower 
proper motion. For hypothetical proper motions 
that are higher, the hypothetical event will be “sped 
up” (tz, ter, and tą will all be shorter) compared 
to the real one, so AxZuicce and Ax, will both fall 


inversely as Moy and so will eventually fall below our 
thresholds of selection and/or detection. On the other 
hand, as Hre] is decreased, the event will be “slowed 
down”, thus reducing thresh in direct proportion, i.e., 
[thresh X fre). Thus, it will eventually fall below our 
search limit of [thresh > 1 masyr™t. Our main concern 
here will be the proper-motion selection effects for the 4 
FFPs. These have 0g = (4.35, 4.77,5.94, 8.80) pas, 


[rel = (4.17, 6.14, 7.95, 9.74) masyr7!, 
Hthresh = (3.46, 11.17, 10.34, 8.80) masyr~!, 
Ax2rect = (2180, 2663, 2992, 6249), and Ax2,) 


(577, 234, 187,1513). Therefore, they have, respec- 
tively, Hrelmin = (1.21,0.55,0.77,1.11)masyr™! and 
remax = (9.09, 16.35, 23.79, 60.87) masyr™!. The 
lower limits play very little role, but the upper limits 
will be important when evaluating the kinematic 
constraints on the FFP masses (see Section). 


8. CONSTRAINTS ON THE FFP POPULATION 


If the “Einstein Desert” is real, then it implies that the 
objects that lie below this gap are a separate popula- 
tion of low-mass objects. The existence of this gap will 
then be an important constraint on the mass function 
of these objects. Hence, the first question that must be 
addressed is, how statistically secure is this gap? 


8.1. Is the Einstein Desert Real? 


The fact that the eye is drawn to the gap in Figure 4l 
does not, in itself, make it real. We can construct 
the following “naive test” by noting that there are 4 
detections within a factor of 2 below the gap, and 6 
within a factor of 2 above it. Given that the sensitiv- 
ity of the survey is monotonically increasing over this 
entire range and that the gap itself covers a factor of 
3, we would expect a continuous distribution to gener- 
ate (4 + 6)/2 x log(3)/log(2) = 7.9 detections in the 
region of the gap. The fact that none are detected has 
a formal probability of p = 0.0004. However, we call 
this test “naive” because it is constructed after the fact. 
The above statistical evaluation would be valid only if 
we asked the question before conducting the experiment. 
We might have asked a dozen questions, such as: is there 
a gap around g ~ 15 pas, around 30 pas etc. Or per- 
haps the data would have yielded some other peculiar 
feature, causing us to construct some other posterior 
test designed to highlight its improbability. If we could 
imagine constructing 50 such tests, then the probability 
of finding p = 0.0004 in one of them is 2%. This is still 
small but the case would not be overwhelming. 
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Nevertheless, the fact is that we did not enter this 
investigation without any prior knowledge. 
already found an analogous gap in the Einstein 
timescale distribution at tz ~ 0.5 days. At typical lens- 
source relative proper motions Hre ~ Tmasyr~! (for 
disk lenses), we would expect a gap at 0p = Hreate ~ 
10 was, which is similar to the geometric center of the 
gap in Figure Æ} 15 pas. Hence, the above statistical 
test should be taken approximately at face value. That 
is, we regard the suggestion by of a 


new population as confirmed}. 


8.2. 6-Function FFP Mass Function 


When and first 
suggested their respective FFP populations, they each 
presented their frequency estimates in terms of ð- 
function mass functions. Neither argued that the FFP 
mass function was sharply peaked, but rather used 6- 
functions as a convenient way to characterize the basic 
mass scale and frequency. In Section we will argue 
that, especially in light of the Einstein radius measure- 
ments presented here, power laws provide a more useful 
framework to characterize FFPs. However, there is one 
constraint on the FFP mass function that is most eas- 
ily described in terms of 6-functions, namely, that the 
observed FFP proper motions put a lower limit on the 
FFP mass scale. 

As noted in the discussion of Figure [5] the observed 
distribution is consistent with that of typical disk lenses. 
Strictly speaking, however, this really applies only to 
disk lenses with D,=22kpc for which the component of 
the proper motion that is due to the peculiar motion of 
the lenses (relative to Galactic rotation) is small com- 
pared to the mean relative proper motion that is set by 
the proper motions of the bulge sources. For Dz, <1 kpc, 
by contrast, these peculiar motions dominate, and the 
expected amplitude of the relative proper motion grows 
Hrel X De with declining distance. Hence, at suf- 
ficiently low M (hence, low Dz), the observed “typi- 
cal disk” proper motions of the FFPs should come into 
strong tension with what is expected for nearby lenses. 

We quantify this conflict as follows. For each as- 
sumed lens mass M, and for each of the four FFP 
events, we calculate Tr. = 62,/KM. We adopt ms from 
and so find Dr = AU/(mel + 7s), 
and we adopt Hs he) from Gaia. We model the disk 
as having a flat rotation curve with vrot = 235kms7!. 
We model the disk lenses as having velocity disper- 
sions o; = 28kms—',/7 and op = 18kms~—',/7, where 
n = exp(D,/2.5kpc), and we assume an asymmetric 


3In fact, (2017) adopted exactly the cautious ap- 
proach recapitulated above toward the evidence that they pre- 
sented for a new population. The statistical significance for 
this suggestion was roughly comparable to ours: they found 6 
events below the gap (compared to our 4). However, with less 
complete light-curve coverage and no finite-source effects, these 
were individually less secure. Their (very appropriate) conserva- 
tive orientation was reflected in their title “No large population 
of unbound or wide-orbit Jupiter-mass planets”, which did not 
mention a new population. The abstract introduces this sug- 
gestion with the cautionary phrase “may indicate”. 


drift of Vrot Us, — (18? + 282 + 332)n (kms—t)?. We 
conduct a 2-dimensional integral over this velocity dis- 
tribution, first calculating H,elne] and than converting 
to We At that point, we exclude realizations that ex- 
ceed the upper limits that are described in Section [7] for 
each event. And, of course, we weight the result by Hrel- 
We find the fraction of this predicted Hre; distribution 
that have proper motions less than the observed values, 
gi. For a fair sample, we expect that the g; should be 
uniformly distributed, with mean (g;) = 1/2, but if the 
model is systematically overestimating the proper mo- 
tions (i.e., it assumes that the lenses are closer than they 
actually are), then we expect g; K 1/2. Hence we adopt 
a likelihood estimator L = i gi. We find that FFPs 
with M < 1.5 Mẹ are disfavored at 2 ø, while those with 
M < 0.95 Me are ruled out at 30. These limits corre- 
spond to mean lens distances of 1.1 kpc and 0.7 kpc, 
respectively. That is, these results conform to our naive 
expectation. 

Thus, any compact mass function (not necessarily 
a 6-function) for which the expected detections were 
MsMe would be heavily disfavored. However, as we 
will see explicitly in Section[8.3] only a small fraction of 
expected detections from viable power-law distributions 
are in this range. 


8.3. Power-Law FFP Mass Function 


The FFP mass function is likely to be better character- 
ized by a power-law function, 


dN M E 
———— = ; (9) 
dlog M Meern 


than a 6-function. Like the latter, a power-law function 
requires two parameters, namely the power-law index 
p and the normalization Z, while the reference mass 
Mporm is an arbitrary zero point whose inclusion allows 
Z to have simple units; i.e., (dex)~!. The first argu- 
ment favoring a power-law distribution is that there is 
no reason to expect that FFPs will have special mass 
scale, only that selection biases may possibly “pick out” 
FFPs of a certain scale. Second, FFPs are very likely 
to be related in some way to bound planets, whose 
mass-ratio function is very well characterized by a power 
law. For lenses with an observed 0p, the inferred mass 
varies as a function of distance, M = 62, /KT rel, with 
Trel = AU(D;" — ie), and so can be very different if 
the lenses are assumed to be in the disk or bulge. Only 
a continuous mass distribution (of which, a power-law 
function has the simplest form) can handle these possi- 
bilities simultaneously. 

Moreover, we do in fact have some prior informa- 
tion on the power-law index p. For planets in the mass 
range under consideration, M < Mgaturn, it is difficult 
to imagine any formation scenario other than within 
the disks around stars, i.e., the same process that gives 
rise to the known population of bound planets found by 
microlensing planet searches. These have a power-law 
mass-ratio distribution, dN/dlogq x q~7, with y ~ 0.6 
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Figure 7. Cumulative distributions of FFP detections for 
three power-law indices (colored curves) compared to the 
actual detections (black). The curves are each normalized 
to 4 detections at Einstein radius 6g = 9 pas. The three 
curves then predict (3,1.5,1) detections within the observed 
Einstein Desert for p = (0.6, 0.9, 1.2), which argues strongly 
that the index is steeper than the measured one for bound 
planets, y = 0.6. 
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Figure 8. Relative number of FFP detections predicted for 
p = 0.9 and p = 1.2 power-law distributions, as a function 
of distance from the Sun and FFP mass. The contours rep- 
resent steps of a factor 1.5. For p = 1.2, the larger area 
in yellow and green almost perfectly compensates the lower 
amplitude compared to the smaller black region, so that the 
FFPs are predicted to come equally from the bulge and the 
disk. For p = 0.9 the ratio is bulge:disk = 5:4. 
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Figure 9. Cumulative distribution of predicted FSPL detec- 
tions (colored curves) versus the actual detections (black) 
of 29 events. The FFP model (p = 1.2), which flattens 
for log(@z/pas)>1.2, is normalized to 4 detections, shown 
as disk (blue), bulge (red), and combined (magenta) curves. 
The same three colors are used for the predictions for stars 
and BDs, which start to rise at log(@z/pas) ~ 1.2. The 
star+BD normalization is set so that the overall prediction 
(green) matches the observed profile. 


(Jung_et_al., in preparation, Zang et _al., in prepara- 
tion, Buzuki et al 2016 Shvartzvald et al][2014). The 
FFPs must either be drawn from this population and 
ejected to larger and/or unbound orbits, or they must 
form in the outer portions of their solar system, i.e., 
beyond the range of detection of the known bound pop- 
ulation. In either case, we would expect this to result in 
a mass distribution with p > y. That is, the probability 
for a planet drawn from the known bound population 
to be scattered to a large or unbound orbit monotoni- 
cally increases with declining mass because this gener- 
ally requires scattering off a heavier planet. And planets 
formed in the outer disk (by core-accretion rather than 
gravitational collapse) should be heavily biased toward 
low mass because of the paucity of raw material and the 
long dynamical timescales. For example, in our own So- 
lar System only dwarf planets such as Pluto are believed 
to have formed in situ in the outer Solar System, while 
Uranus and Neptune are believed to have formed inside 
the current orbit of Saturn and to have been scattered 
to their present-day orbits. 

In the observed mass-ratio function for bound plan- 
ets, the number of planets with mass ratio greater than 
a certain q, is x q`. Then, assuming that scattering 
to large or unbound orbits is proportional to the num- 
ber of such potential scatterers, p = ky, with k = 2. 
In practice, heavier planets will be more efficient scat- 
terers, which would imply 1 < k < 2. On the other 
hand, if the number of FFPs is of order or larger than 
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the number of bound planets (as seems to be the case, 
see Section [9.2), then the original value of y must have 
been larger than currently observed, with lower mass 
members of the original bound population having been 
preferentially removed. In this case, the present-day ob- 
served mass ratio distribution (and, in particular, its 
power-law form) must reflect the imprint of the scatter- 
ing process in addition to the formation process. How- 
ever, because the FFP population is not overwhelm- 
ingly dominant, it is unlikely that taking account of 
formation-plus-scattering can drive k much above the 
rangel<k <2. 

Starting from these very general considerations, we 
ask whether and how the FSPL FFP sample can fur- 
ther constrain p. Figure [Z]shows the predicted cumula- 
tive distribution of FFPs based on the range of power 
laws described above, p = (0.6,0.9,1.2), correspond- 
ing to k = (1,1.5,2). In making these predictions we 
have adopted relative detection efficiencies € = (4/3)X 
(X < 0.5), E = (2X + 1)/3 (05 < X < 1), E= 1 
(<X<2), € =2-xX (2 < X < 3), with 
X = log(fg/3 pas). In the regime X > 0.5, in which 
the detections are roughly randomly distributed in the 
range of detectability (Figure [6), this function approx- 
imates the relative fraction of events that are sensitive 
to a given 0g. For X < 0.5, we complete this function 
linearly by imposing a threshold at 6g = 3 uas, which is 
supported by the fact that all four FFPs are pressed up 
close to this limit. 

The curves are all normalized to the 4 detections with 
0g < 9 uas. As can be seen, all three curves account of 
the form of the detections very well. However, they then 
predict, respectively, about 3, 1.5, and 1 event in the 
Einstein Desert. Hence, they are, respectively, strongly 
disfavored, mildly disfavored, and acceptable. 

In Figure[8] we show the predicted distribution of de- 
tections as a function of FFP mass and distance for 
the latter two models. The contours change color at 
factors of 1.5. The bulge FFPs are centered at M ~ 
(40,30) Me for p = (0.9,1.2), while the disk FFPs are 
at a broad range of lower masses, as expected from the 
facts that 0g is measured for the objects in our sample 
and M = 02, /KTre1 As we will soon see, there are al- 
most exactly the same number of bulge and disk FFPs 
in the p = 1.2 model: the fact that the yellow and green 
contours are much larger than the black contour com- 
pensates for their 2.2-3.4 times lower amplitude. For 
the p = 0.9 model, we find that the ratio is bulge:disk 
= 5:4. 

To normalize this distribution we predict the stellar 
and brown-dwarf FSPL events using a Chabrier initial 
mass function (with bulge stars that have M > 1Mo 
being converted into remnants). See Figure[9] We nor- 
malize the FFP curves to 4 total detections and nor- 
malize the stellar curve to approximately match the ob- 
served distribution. The bulge and disk FSPL events 
are shown in red and blue respectively, with magenta 
showing the sum of these. There are two sets of such 
curves, one for FFPs and the other for stars and BDs. 
The green curve shows the overall sum. For the FFPs, 


the bulge and disk curves are almost identical. As noted 
above, for the p = 0.9 model, the bulge FFPs would be 
favored 5:4. From the relative normalization of the FFP 
and star-plus-BD curves, the FFP mass function can be 
expressed as 


dN 0.39 ( M 


=p 
ee = — = 0.9 or 1.2). 
dlogM dexx star \ 38 <) » orka] 

(10) 
Because the mean mass of the adopted stellar-BD mass 


function is (M) = 0.199 Mo, this can also be written 


dN 1.96 M P 
—— = —— x | —— ; (11) 
dlog M dex x Mọ stars 38 Me 


The uncertainty in the normalization is strongly domi- 
nated by the Poisson error of 4 FFP detections, i.e., it 
is 50%. We consider that it is unlikely that the power- 
law index lies far outside the indicated range 0.9Sp<1.2. 
For indices much below p = 0.9, the FFPs would pop- 
ulate the Einstein Desert. While there are no strong 
constraints on power laws that are substantially steeper 
than p = 1.2, it is hard to imagine a mechanism that 
would generate them. 

Assuming that these power laws apply to FFPs from 
Earth to Saturn masses, then this range contains a total 
of 130 Mẹ or 200 Mẹ of FFPs per solar mass of stars for 
the p = 0.9 and p = 1.2 models, respectively. 

We note that the form of the bulge and disk stel- 
lar FFP detections shown in Figure [9] are in reason- 
able accord with the expectations of Section B] The 
observed FSPL events are above the expectations in the 
BD regime. This may reflect an underestimation of the 
BD component in the Chabrier mass function, but it 
could also reflect a statistical fluctuation. Addressing 
this question would take us beyond the scope of the 
present investigation. 


9. DISCUSSION: RELATION OF FFPS TO KNOWN 
POPULATIONS 


When we introduced the acronym “FFP”, we were care- 
ful to have it refer to “candidates”, rather than just 
“free-floating planets”. The main reason for this is the 
general concern that the low-mass objects that appear 
“isolated” in microlensing events could be bound planets 
in orbits that are too wide to enable microlensing signa- 
tures of the host. This possibility has always been rec- 
ognized and was explored in depth by (Clanton & Gandi 
Poal, (2016). paed aai hae thetor 
bital separations of FFPs in very wide (“Kuiper-like” ) 
and extremely wide (“Oort-like”) orbits could eventu- 
ally be measured. 

Here, we compare the FFP population as derived from 
our FSPL survey, as characterized by Equations (0) 
and (lI) to various known populations, in order to ob- 
tain a more comprehensive picture. This will include 
planets discovered by microlensing in typical orbits, 
known planets in very wide orbits, and known inter- 
stellar objects. 
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We adopt the orientation that the most likely origin of 
FFPs is that they they formed within a factor of a few of 
the snow line, a region that was rich in proto-planetary 
materials, and where dynamical timescales were rela- 
tively short, and that they were ejected by dynamical 
processes, either to much wider orbits or to unbound 
orbits. Hence, we regard it as likely that FF Ps comprise 
both bound and unbound objects. Within this frame- 
work, we expect a greater fraction of high-mass FFPs to 
be in bound orbits, simply because it is easier to eject a 
test particle than an object whose mass is a substantial 
fraction of that of the perturber. 


9.1. Comparison to Mroz et al. (2017) 


considered a 6-function model with 
FFP mass M = 5 Mg to account for their Nip = 6 short- 
tg events. They obtained a best fit of Rs Mọ,tp = 10 
FFPs per star. Although they did not quote an er- 
ror in this estimate, the Poisson noise yields o5 mg ,ty = 


Rs Me te / NEP = 4.1. 

To compare our result to theirs, we insert the same 
-function model into the formalism described in Sec- 
tion[8.3] thereby obtaining Rs wz,o, = 15.3 7.6, based 
on our Nọ, = 4 small-ĝg detections. Hence, the two 
determinations are consistent at the 0.6 ø level. 

Although we mainly restrict consideration to our own 
results in this work, we note that, because the two nor- 
malizations are consistent, it can also be appropriate to 
combine them. We then find that the combined normal- 
ization is a factor 


Frenorm = 0.73 jæ anj 0.24 


(12) 


smaller compared to results from our study alone. Note 
that the fractional error is then ~ 1/3 rather than 1/2, 
as would be expected from the fact that there are a 
total of 10 detections, rather than 4. Note also that this 
renormalization would apply to both Equations (10) and 
and to either indicated power-law index. 


9.2. Comparison to the Known Bound-planet 
Mass-ratio Function 


At present, the masses of most microlensing planets are 
unknown. Hence, what is most precisely measured is 
the planet-host mass-ratio (q) function, which can be 
expressed, in analogy to Equation (9), 


dN a 
dlogq a x ( ) i 
d log q norm 


which, using a zero point qnorm = 10735, has measured 
values Z, = 0.175 (dex)~?/star and y = 0.6 (Jung et 
al. 2022, in preparation, Zang et al. 2022, in prepara- 
tion). There are three things that make the comparison 
difficult. First, the units of Z are (dex)~'/star while 
the units of Z, are (dex)~?/star. The additional “dex” 
means “per decade of projected separation”. We ac- 
count for this by regarding two decades of separation 
(roughly 0.1-10 AU) as being approximately represen- 
tative of the whole population. That is, we roughly 


(13) 


account for not only the cold planets detected by mi- 
crolensing but also the hot and warm planets detected 
by transit and radial velocity surveys. Second, the de- 
nominator term “star” does not mean exactly the same 
thing for the two power-law functions. For our FSPL 
study, events enter the sample in proportion to their 
number density, but for the mass-ratio functions studies 
of Jung et al. (in preparation) and Zang et al. (in prepa- 
ration), they enter as microlensing events, whose fre- 
quency is weighted by their cross section, i.e., x M1/?. 
Therefore, if we want to convert Equation (13) to “per 
solar mass of stars”, we should calculate the mean mass 
weighted by M1/?, i.e., (M) q7 ~ 0.33Mo. Third, 
eponymously, the argument of the mass function is mass, 
while the argument of the mass-ratio function is mass 
ratio. We account for this by evaluating at bound-planet 
mass Mp = q(M) qq. With these approximations, we 
infer a frequency of known bound planets of 


aN 1 (me Ba (14) 
dlog Mp 7 dex x Mo,stars 38 Me i 


Comparison of Equations and indicates that 
at the normalization mass there are 2 times more FFPs 
than bound planets, while at lower masses, the ratio is 
even larger. For example, at mp = 10 Mg, it is about 
3 for the p = 0.9 power law. Of course, this compari- 
son required several approximations due to the intrin- 
sic incommensurablity of the underlying measurements. 
Nevertheless, one can infer robustly that the number of 
FFPs is of order or larger than the number of bound 
planets in the mass range to which the FFP measure- 
ment is directly sensitive, 5SM/Me<60 (see FigureB). 
This implies, as foreshadowed in Section [8.3] that the 
form of the bound-planet mass function was shaped 
partly, or perhaps mainly, by the ejection process, rather 
than just the formation process. 


9.3. Comparison to Known Uranus- and Neptune-like 
Planets 


In the entire Universe, there are exactly three known 
planets that can be securely described as “Uranus-like” 
( or “Neptune-like”), in that their masses and separa- 
tions are of the same order as these two planets, namely, 
Uranus and Neptune themselves, as well as OGLE-2008- 
BLG-092Lb, which has s = 5.3 and q = 2.4 x 1074 
(Poleski et all, 2014. 


9.3.1. Uranus and Neptune 


We focus first on Uranus and Neptune. Equation 
predicts that, for p = 0.9 and within a factor v10 of 
M = (Muranus + Mneptune)/2 = 16 Mg, each solar mass 
of stars should be accompanied by an average 4.3 FFPs 
(whether unbound or in wide orbits). This can be com- 
pared to the two such observed planets in the Solar Sys- 
tem, the only planetary system for which we have a com- 
plete census. These numbers are compatible within the 
Poisson error, and they would be more so if we applied 
the renormalization given by Equation (12), by which 
the expectation would be reduced to 3.1 + 1.0. 
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Furthermore, if a doppelganger of the Solar System 
gave rise to a microlensing event at Dz ~ 5kpc in which 
either its “Uranus” or its “Neptune” generated a char- 
acteristic FSPL profile, there is only a small chance that 
its “Sun” would leave a trace on the event. That is, if a 
planet lies at (3-dimensional) physical separation from 
its host, r, then the probability that its normalized pro- 
jected separation, s, and host-source impact parameter, 
ug, lie within specified regimes are, 


2 


S U 
p(s < Smax) =1-,/1- “2 pluo < U0,max) = i , 
Sfull Sfull 
(15) 


where sf = r/Dr0g. We note that for Dg = 5kpc 
and Ds = 8kpc, Dz0g = 3.91 AU(M/Mọo)!/?, but for 
the entire range, 2.5 kpe < Dz < 5.5 kpc, this formula 
remains accurate to < 4%. Thus, for Uranus and Nep- 
tune doppelgangers, sf, = 4.9 and 7.7, respectively. 
We focus first on Neptune. According to Equation 
there is a small (p = 8%) chance that s < 3, in which 
case, even if the host were detected, the planet would 
be regarded as having a Saturn-like (not Neptune-like) 
orbit. For FSPL events, the main way that the host 
would be detected is that the source would come within 
uo < Uomax ~ 2.5, which would induce a peak mag- 
nification A21.03. Overall, this occurs with probabil- 
ity p = 32%, but about 7% is due to the cases with 
s < 3, leaving p = 25% that could in principle be rec- 
ognized as having Neptune-like orbits. At the median 
projected separation of the s > 3 population, namely 
Smed = 6.8, the encounter with the host would occur 
At ~ (Smed/ Stun)? / Dr bre — 68 day (frei /6 mas ce a 
before or after to. To reliably identify a peak in such 
a low-amplitude event, we should require that the ob- 
servations continue a factor 1.3 beyond peak, i.e., to 77 
days. Roughly 1/3 of all events will fail this criterion be- 
cause the peak region lies wholly or partly outside the 
microlensing season. Thus, we expect about 16% of true 
bound Neptunes that are detected as FSPL FFP events 
to be recognizable as such. Repeating the same calcu- 
lation for Uranus yields 32%. Thus, the average prob- 
ability for Uranus and Neptune is about 25%. Hence, 
with four examples, we should expect one host detection 
assuming that all were bound. However, the probability 
for zero detections under this assumption is p = 32%. 
Moreover, there is added uncertainty in this estimate 
due to the fact that the majority of potential hosts are 
much less massive than the Sun. If typical separations 
of Neptune-like planets scale x M1/?, then the above 
calculations hold exactly. If they they scale x M, then 
the probability of host detections goes up, while if these 
typical separations are independent of host mass, then 
the probability goes down. 


The main implication of this exercise is that, at least 
regarding the region of parameter space in which Equa- 
tion is directly sensitive to the observations, the 
Solar System is consistent with the hypothesis that all 
FFPs are due to bound planets in wide orbits, such as 
those of Uranus and Neptune. 


9.3.2. OGLE-2008-BLG-092Lb 


Among bound planets with secure s and q measure- 
ments: OGLE-2008-BLG-092Lb is unique: there are no 
other planets with s > 3 and q < 10718. Planets with 
s < 3 are unlikely to be in wide orbits, while those with 
q > 107Ł8 could in principle have formed through an- 
other channel, i.e., gravitational collapse. 

However, in addition to this one secure Neptune- 
like planet, there is another bound planet, OGLE-2011- 
BLG-0173Lb, with best estimates s = 4.7 and q = 
4.5 x 10-4 (2018). In fact, there is an 
alternate solution with s = 0.22 with virtually identi- 
cal y?. Nevertheless, based on the general statistical 
properties of microlensing planets, (2018 
argued that the wide solution was preferred. 

conducted a systematic search for 
wide-orbit planets in 20 years of OGLE data, which 
recovered these two planets (and no others in the 
“Neptune-like” domain). They re-evaluated and con- 
firmed the argument that the wide solution for OGLE- 
2011-BLG-0173Lb is preferred, based on updated mi- 
crolensing statistics. Including both planets, they de- 
rived a rate of 1.41} “wide-orbit ice giants” (2 < s < 6, 
—4 < logg < —1.5), based on a total of 5 detections, 
only the above two of which are in the regime of inter- 
est here. 

To obtain a proper comparison to the FSPL sam- 
ple, we apply their best-fit parameters (A,n,m) = 
(1.04, 1.15, 1.09) to their Equations (6) and (8), but re- 
stricted to (4 < s < 6, —4 < logq < —3), and find an 
expected number (per star) of N = 0.67+0.47, where we 
have adopted a Poisson error based on two detections. 
To maintain consistency with our other comparisons, 
we should convert this to expected Neptune-like plan- 
ets per solar mass of stars, i.e., 2.0 + 1.4. The best-fit 
value is identical to the actual number of such planets in 
the Solar System, and, additionally, there is a substan- 
tial Poisson error in this prediction. Hence, this result 
is likewise consistent with all Neptune-like FFPs being 
due to wide-orbit bound planets. 


9.3.3. Caveats 


Although both of the above tests indicate consistency 
with the hypothesis that all of the FFPs in our FSPL 
sample are due to known populations of wide-orbit 
bound planets, there is no strong evidence that even 
a majority (let alone all) are due to such planets. First, 
the solar system is known to be rich in gas giant planets 
in “Jupiter-like” or “Saturn-like” orbits relative to the 
field. From the mass-ratio function, we would expect 
only 0.2 such planets within v10 of their mean mass, 
compared to two observed. Moreover, the presence of 
Jupiter and Saturn is significant not only in that they 
may indicate that the Solar System is planet-rich, but 
also in that they were implicated in the expulsion of 
Uranus and Neptune to wide orbits. Finally, the ob- 
served characteristics of the Solar System may be sub- 
ject to subtle selection effects due to the existence of the 
observers. For example, the turmoil in the Solar System 
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engendered by the expulsion process may have played a 
crucial role in the delivery to Earth of water and/or or- 
ganic materials, or perhaps to other conditions on Earth 
that were crucial to the development of intelligent life. 

Similarly, the estimates of Neptune-like planets from 
the detection of wide-orbit planets are subject to sev- 
eral caveats. First, the Poisson error in this estimate is 
already quite large. Second, one cannot really be cer- 
tain that the wide-orbit solution for OGLE-2011-BLG- 
0173Lb is correct. This solution is favored statistically 
only because there are more planets with s > 1 than 
s < 1, but the power law relation is derived from a 
sample that does not contain any planets (other than 
OGLE-2008-BLG-092) that are close to its observed 
value: |logs| ~ 0.7 and that are also logq < —1.8 (so 
not subject to alternative formation mechanisms). If 
this planet is eliminated, then the mean expected num- 
ber drops by a factor 2 and the fractional Poisson error 
becomes much larger. 

What we can say robustly is that there is a known 
population of bound planets that contributes signifi- 
cantly to FSPL events that have been (or possibly will 
be) in the FFP domain. However, whether this is 10% 
or 100%, we cannot now say. 

Finally, we note that (contrary to the case of the Jo- 
vian FFPs suggested eta there are no 
limits on bound Neptune-like planets from direct imag- 
ing because they are too faint. 


9.4. Comparison to ?Oumuamua-like Objects 


’Oumuamua (Williams, was discovered via the 
PanStarrs survey as an asteroid-like object that was 
leaving the Solar System in a strongly unbound orbit 
(eccentricity e = 1.2). It may or may not be related 
to FFPs, depending on its physical origin. The conven- 
tional view, and the one that we will adopt here, is that 
it is an interstellar asteroid that was formed in another 
solar system. However, other ideas for its origin have 
been suggested, including that it is an alien spacecraft 
(Loeb, and that it is the result of a violent col- 
lision of Solar System objects that occurred inside the 
orbit of Mercury (B. Katz et al, 2018). 

calculated the space density of ob- 
jects with exactly the same physical properties as 
’Oumuamua under the assumption that the expected 
rate of detection was exactly 1 per 3.5 years of PanStarrs 
operations as n ~ 0.20 AU~® = 1.75 x 107 pe. To 
account for the fact that there have been no further re- 
ports of such objects during the past 4 years, we reduce 
this to n = 1 x 10715 pc. Following [Rafikov 
we adopt physical dimensions (180mx18mx18m), and 
a typical asteroid density of 3 gm/cm? to derive a mass 
Moum = 1.7 x 10". gm = 3x 1071!" Mg, and so a spatial 
density Poum = 0.03 Mg pe~?. In order to compare this 
to our FFP estimates, we first adopt a local stellar mass 
density of pstars = 0.03 Mopc™? and assume that this 
detectability applies to 1 dex of object mass, to obtain 
dNoum — 3.3 x 1016 , (16) 
dlog M dex x Mo stars 


implying an asteroid mass density per decade of 
1.0 Me/dex/Mo stars: 

The main point to notice about this calculation is 
that this mass fraction per decade (measured at M = 
3 x 10717 Mo) is remarkably similar to the value im- 
plied by Equation (II) at its pivot point (M = 38 Mẹ) 
of 32 Ma/dex/Mo stars. That is, the two are consistent 
with being part of the same distribution provided that 
the distribution function is approximately flat in mass, 
i.e., p~ 1. As we have already concluded that p ~ 1 
based on FFP detections, combined with the lack of de- 
tections in the Einstein Desert, this suggests that FFPs 
and ’Oumuamua-like objects may result from the same 
physical processes, i.e., the processes of planet formation 
and early dynamical evolution. 

In fact, considered and rejected the 
idea that "Oumuamua-like objects were created and 
ejected as part of planet formation, primarily because 
they appeared to be too abundant. However, we con- 
sider that the error in the above estimate is an order of 
magnitude (i.e., 1 dex). In particular, if ?Oumuamua- 
like objects were a factor 10 less common than the 
estimate, there would still be a 10% 
chance that one would be detected, which hardly can be 
considered to be an implausible scenario. 

If we assume that they are part of the same process 
and that they participate in the same power-law distri- 
bution, then by comparing the two measurements sepa- 
rated by 18.1 dex in mass, which differ in amplitude by 
—1.5 + 1 dex, we can estimate a power-law index of 
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However, it is also possible that the similarities of the 
mass densities of FFPs and ’Oumuamua-like objects is 
just a coincidence and that ?Oumuamua arose from some 
entirely different process. In addition to the possibili- 
ties mentioned above, [Rafikov] has suggested that 
*Oumuamua results from the tidal disruption of a rocky 
planet by its white-dwarf host. There could be other 
possibilities as well. 

Nevertheless, the hypothesis of a common origin can 
be subjected to further tests. If more ’"Oumuamua-like 
objects are found by improved surveys, their frequency 
and power-law distribution can be better measured. If 
this power law is consistent with the p ~ 1 that is im- 
plied by the common-origin hypothesis, this would con- 
stitute further evidence. Similarly, better measurements 
of the FFP power law could confirm (or contradict) its 
apparent agreement with the one needed to connect the 
FFP and ’Oumuamua measured frequencies. 

In this regard, it is important to note that there 
has already been a detection of a second extra-solar 


minor body, 2I/Borisov (Higuchi & Kokubo, (2020), of 
substantially different mass, M ~ 5 x 10714 Mẹ 
(Jewitt et all 2020, i.e., 1700 times more massive than 
‘Oumuamua. Unfortunately, to our knowledge, there is 


not as yet a published estimate of the space density of 
objects of this class. 
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9.5. Summary 


The comparison to various known populations can be 
summarized as follows. First, the frequency of FFPs de- 
rived from our FSPL sample is consistent with that from 
the PSPL sample of [Mréz et all at the 1ø level. 
Combining the two leads to pre-factors of 0.28 + 0.09 
and 1.43 +0.47 in Equations (10) and (Ii), respectively. 
Second, the frequency of FF Ps is of order or larger than 
that of known bound planets. Third, if the Solar Sys- 
tem’s endowment of planets with similar masses and or- 
bits to those of Uranus and Neptune, then these objects 
are likely to contribute a substantial fraction of our de- 
tections. However, there is some evidence that the Solar 
System may not be typical. Fourth, the observed FFPs 
and ’Oumuamua-like objects are consistent with being 
drawn from the same power-law distribution, in which 
case its index would be p = 0.92 + 0.06. 


APPENDIX A. REMARKS ON INDIVIDUAL EVENTS 


In this section, we remark on anything that is notable 
about the FSPL events from 2016-2018. For notes on 
the 2019 events, see [Kim et al] (2021). 

The majority of these notes concern the fact that 10 
of the 17 FSPL events that are reported here were previ- 
ously published or (in one case) is the subject of work in 
preparat trast, only one of the 13 events ana- 


ion. By con 
lyzed by (Kim eral from 2019 had been the sub- 
ject of a previous publication, namely the FFP OGLE- 
2019-BLG-0551 (Mróz et all [2020). For each of the 9 
published events, we compare the best value of 0g re- 
ported in these papers (first) to the one we report here 
(second), both in pas: 

Two of these events, OGLE-2016-BLG-1540 
(Mróz et_alJ,|2018) (9.2 vs. 8.8) and KMT-2017-BLG- 
2820 (Ryu et al], (5.9 vs. 5.9) are FFPs, with the 
latter having been discovered as part of the program 
being reported here. 

Three were published as BD candidates as inferred 
from their relatively small Einstein radii: OGLE-2017- 
BLG-0560 (Mróz et al] 39 vs. 35), MOA- 
2017-BLG-147 (Han et al, (2020) (51 vs. 46), and 
MOA-2017-BLG-241 (Han et al.) |2020) (28 vs. 26). 

Four of these events were the subject of previous pub- 
lications because they are FSPL events for which it was 
possible to derive microlens parallaxes based on Spitzer 
data. Spitzer carried out a large, 6-year program whose 
principal goal was to measure the microlens parallaxes, 
and thereby the masses of and distances to, planetary- 
system microlenses (Yee et al], 2015). However, Spitzer 
microlens parallaxes also yield masses and distances for 
FSPL events, whether the finite source effects are ob- 
served from the ground or Spitzer (Zhu et al], (2016). 
The four published events are: OGLE-2016-BLG- 


1045 (Shin et al} |2018) (244 vs. 249), OGLE-2017- 
BLG-0896 (Shvartzvald et al), (2019) (140 vs. 140), 
OGLE-2017-BLG-1186 (Li et al], 2019) (94 vs. 94), 
and OGLE-2017-BLG-1254 ERTS Sa (207 


vs. 224). 
The reason that the 0g values are identical for three 


events (KMT-2017-BLG-2820, OGLE-2017-BLG-0896, 
OGLE-2017-BLG-1186) is that we adopted the pub- 
lished parameters after finding that we recovered the 
event using our standard procedures. 

In addition, S. Tsirulik, J.C. Yee et al., (in prepa- 
ration) applied this technique to OGLE-2018-BLG- 
0705 and found Dz = 2.2 + 0.1 kpc and M = 0.74 + 
0.05 Mo. The event is projected against the cluster 
NGC 6544, and the proper motion measured by Gaia 
of the blended light (almost certainly the lens) is consis- 
tent with that of the cluster. Hence, it is very likely that 
the lens is a cluster member, having 2 ø tension with the 
literature-based distance to NGC 6544 (2.6 + 0.3 kpc). 

OGLE-2017-BLG-0905: This event is by far the 
largest outlier on the CMD, lying roughly 2 mag above 
the bulge giant branch (highest point in Figure). This 
“problem” would be ameliorated (but not completely 
solved) if the true color were substantially redder than 
the one that we measured. First, however, we have 
made separate color measurements from KMTC and 
KMTS data, and these differ by only 0.04 mag. Sec- 
ond, the color (and the magnitude) that are derived 
from the light curve are nearly identical to those of the 
baseline object, as would be expected for such a bright 
source. This is true for both KMTC and KMTS. An- 
other possible explanation for the extreme brightness of 
the source (for its color) is that it is actually a fore- 
ground disk star. However, the Gaia proper motion 
us(E, N) = (—5.3 + 0.3,-7.5 + 0.2)masyr7!, is only 
~ 3 mas yr™! from the centroid of the bulge distribution. 
On the other hand, this zg would be unusual for a disk 
star, unless it were very nearby, which is disfavored by 
the Gaia parallax measurement 7g = 0.13 + 0.27 mas. 
Thus, the source is most likely a relatively rare post- 
AGB star. 

Although ug is only lo from the bulge mean, this 
deviation happens to be close to the direction of anti- 
rotation, so that us = 9.2masyr~'. This makes 
our otherwise slightly puzzling measurement of [ye = 
12.1 masyr™! (the second largest in our sample), more 
understandable. 

Finally, we note that most of the J-band points in and 
near the peak of the event were saturated. However, 
the event parameters were easily recovered from the V- 
band light curve. Moreover, a significant minority of I- 
band points were not saturated due to a combination of 
“poor” seeing and/or favorable placement of the source 
near the pixel corners. These “salvaged” J-band points 
yielded consistent results with the V-band analysis. 

Two other events, MOA-2016-BLG-258 and 
OGLE-2017-BLG-1186, were also saturated at peak. 
In the first case, the saturation was mild and intermit- 
tent, and saturation did not affect KMTA at all. Hence, 
no special measures were required to fit the light curve 
other than removing a few saturated points. For the 
second case, which was already mentioned above as a 
published event, the peak region was very densely cov- 
ered in the V band due to a combination of the long 
source self-crossing time tą = pte = 3.8days and the 
high V-band cadence I = 0.4 hr’. Hence, there was no 
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difficulty in modeling the event even though all of the 
peak J-band data had to be excluded. 

OGLE-2017-BLG-0084: This event has a relatively 
large color error, o(V — I)g = 0.14, in a measure- 
ment that yields a somewhat unexpectedly red source, 
(V —I)so = 1.33 + 0.14, compared to (V — I)so0 ~ 
1.06 + 0.10 for lower red-giant branch stars of its appar- 
ent magnitude and in its angular proximity. The rela- 
tively poor measurement is due to the event peaking at 
the beginning of the season, when the nightly visibility 
was brief, so few V-band points were taken when the 
Vs ~ 23.5 source was highly magnified. This color is 
very similar to that of the baseline object, whose error 
is also large. We accept the color measurement at face 
value. Nevertheless, we note that if we had imposed 
the lower giant-branch color, then the estimates of 0g 
and ure] would have been reduced by a factor 0.83 from 
264 to 211 mas, and from 1.71 to 1.42 mas yr~!, respec- 
tively. We do not impose this prior because it would 
be inappropriate to do so without also taking account 
of the reduced probability of detection, which in this 
loWw-Hre] regime scales roughly as Opy/3.,, i.e., implying a 
combined factor of 0.47. The net effect of applying both 
priors would be too small to warrant introducing this 
level of complexity. 

Three events met the initial selection criteria and also 
yielded excellent p measurements, but were nevertheless 
excluded from the final sample because the sources ei- 
ther are not giants or could not be properly character- 
ized. As a matter of due diligence, we document these 
decisions. 

KMT-2017-BLG-1725: The high magnification of 
this event enables a very precise measurement of the 
(V — I) color, despite high extinction, A; = 3.84, reveal- 
ing that the source lies A[(V — I), I] = (—1.10, +1.03) 
blueward and fainter than the clump. These offsets are 
inconsistent with a bulge giant, and indeed with any 
bulge source, other than from very rare populations, 
such as extreme blue horizontal branch stars. Most 
likely, the source lies far in the foreground along this 
b = —1.1 line of sight, which implies that it would be 
probing a very different population of lenses compared 
to most of our sources. In addition, due to the source’s 
unknown extinction, we cannot reliably estimate 0, (so 
0g) which is the basic point of the entire project. We 
therefore exclude this event from the sample. 

KMT-2018-BLG-1993: This event lies far out on 
the near side of the bar and very close to the Galac- 
tic plane (1, b) = (+6.06, +0.33). The KMT event page 
lists the extinction as A; = 5.60 (derived by approxi- 
mating Ay = 7 Ag, where Ax is ee cornea 
(2019). The KMT V/I CMD confirms that the field 
is too extincted to determine the source properties in 
these bands. As we usually do in such cases, we try to 
clarify the nature of the source by incorporating VVV 
survey data (Minniti et al], (2017). We match KMT I- 
band pyDIA photometry to K-band VVV to construct 
an [J — K),I| CMD. We find that the clump is at 
Ia pypra = 19.90, which combined with Tao = 14.27 
from [Nataf et al] (2013), would yield A; ~ 5.6 (taking 


account of the fact that pyDIA typically differs from 
standard Cousins J by $0.2 mag whenever it can be 
calibrated). From the light curve model, we obtain 
Is pypra = 21.79 + 0.10. implying that the source lies 
1.89 + 0.10 below the clump. In principle, this would 
be consistent with it being a lower giant-branch star in 
the bar. We then check the nearest VVV entry and find 
that it is A[(J— K), K] = (—0.40, 1.51) bluer and fainter 
than the clump. This “catalog entry” cannot be directly 
identified with the microlensed source (with AJ = 1.89) 
because if it were, then it would be A(I — K) = +0.48 
redder than the clump. Therefore, the VVV entry must 
contain a blend that contributes substantial light that 
is much bluer than the source. In principle, the source 
might then still be a giant at the base of the giant 
branch, and, for example, the blend could be the lens, 
lying far in the foreground and therefore blue. Unfor- 
tunately, we have no corroborating evidence that this is 
the case, and the source could be a subgiant lying much 
closer along the line of sight. The plausibility of this 
alternate scenario is augmented by the source’s near- 
plane position, b = +0.33. That is, the line of sight is 
never more than 50 pc from the plane, so without color- 
magnitude information, the source could be anywhere 
along this line of sight. Hence, we exclude this event. 


KMT-2016-BLG-2608: Guided in part by experi- 
ence with the above two events, we carefully reviewed all 
the FSPL events, including those from 2019 that were 
analyzed by[Kim et all (200), in order to identify prob- 
lematic events. This investigation led to the exclusion 
of one further FSPL event: KMT-2016-BLG-2608. This 
event was the 20th “new event” found in our special 
giant-source search, beyond the 2588 events found by 
EventFinder for 2016. Following the conventions es- 
tablished by and (2021), 
we assign it a discovery sequence number 2588 + 20 = 
2608. Due to very heavy extinction, Ay = 5.56, the 
(V — I) color could not be measured despite relatively 
high magnification Amax ~ 40. Like KMT-2018-BLG- 
1993, KMT-2016-BLG-2608 lies very close to the plane, 
b = —0.27, but in contrast to that event, the source 
is much brighter than the clump, i.e., by AJ = —1.38 
mag. In itself, this makes it a very plausible giant-source 
candidate. However, we are unable to obtain mini- 
mal confirmation of this assessment from VVV data. 
The nearest catalog entry has K = 12.58 (compared to 
Ka = 13.80). If we identify this entry with the source, 
then it is A(J— K) ~ —0.16 bluer than the clump. This 
would be slightly unusual for such a bright star, and 
we note in particular that this is not a problem that 
could be resolved by assuming that the VVV entry con- 
tains blended light: this would only make the source 
bluer. These concerns are amplified by the VVV IR 
color information: there is only one other band with a 
flux measurement, H, and this yields (H—K)s = —0.04, 
compared to (H — K) ~ 0.8 for typical bulge giants in 
this field. Because (H — K)o is strictly positive for all 
normal stars, the only plausible explanation for this IR 
color measurement (other than measurement error) is 
that the source lies well in the foreground. This expla- 
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nation is also consistent with the very low Galactic lati- 
tude of this line of sight, b = —0.27. In principle, as just 
mentioned, it is possible that the H-band measurement 
is incorrect, but the giant-source scenario requires that 
this error be about 0.8 mag. While possible in principle, 
the balance of evidence favors a foreground source, so 
we eliminate this PSPL event. 

Finally, we note that according to the above- 
mentioned naming conventions, the KMT counterparts 
to OGLE-2016-BLG-0245 and OGLE-2017-BLG-0560, 
which were both re-discovered in the special giant-source 
searches, are assigned KMT names KMT-2016-BLG- 
2627 and KMT-2017-BLG-2830, respectively. In con- 
trast to KMT-2016-BLG-2608, these names have no 
lasting significance because the official (i-e., discovery) 
name is from OGLE. However, these names are needed 
in the context of the current paper to maintain homo- 
geneous conventions. 

Nevertheless, this purely formal naming issue does re- 
late to an important, substantive question: within the 
sample of 30 FSPL events, exactly four were not discov- 
ered by EventFinder: OGLE-2019-BLG-0551 (p = 4.5), 
OGLE-2016-BLG-1540 (p = 1.6), KMT-2017-BLG-2820 
(p = 1.2), and OGLE-2017-BLG-0560 (p = 0.9). All 
four were discovered in the special giant-source search, 
but the first has a “normal” KMT name because it had 
already been discovered by AlertFinder. This sample 
of EventFinder “failures” all have high p. Indeed, of 
the 5 FSPL events with p > 0.6, only KMT-2019-BLG- 
2073 (p = 1.2) was discovered by EventFinder. Recall 
that the original motivation for the special giant-source 
search was that EventFinder had failed to discover 
OGLE-2019-BLG-0551, likely (it was thought) due to its 
deviations from the standard form of 
and profiles. By contrast, AlertFinder 
simply looks for rising light curves, but this system 
was not operating (or not fully operating) in 2016-2018. 
Hence, the discovery of KMT-2017-BLG-2820 and the 
tight overlap between EventFinder “failures” and high- 
p events together show that this innovation was both 
necessary and (at least mostly) sufficient. 
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